Experimental observation of the optical spin-orbit torque 
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Spin polarized carriers electrically injected into a magnet from an external 
polarizer can exert a spin transfer torque (STT)^on the magnetization. The phe- 
nomenon belongs to the area of spintronics research focusing on manipulating 
magnetic moments by electric fields and is the basis of the emerging technologies 
for scalable magnetoresistive random access memories. 2 In our previous work we 
have reported experimental observation 3 of the optical counterpart of STT^EI in 
which a circularly polarized pump laser pulse acts as the external polarizer, 
allowing to study and utilize the phenomenon on several orders of magnitude 
shorter timescales than in the electric current induced STT. Recently it has 
been theoretically proposed^! an( j experimentally demonstrated^^ that in the 
absence of an external polarizer, carriers in a magnet under applied electric field 
can develop a non-equilibrium spin polarization due to the relativistic spin-orbit 
coupling, resulting in a current induced spin-orbit torque (SOT) acting on the 
magnetization. In this paper we report the observation of the optical counter- 
part of SOT. At picosecond time-scales, we detect excitations of magnetization 
of a ferromagnetic semiconductor (Ga,Mn)As which are independent of the po- 
larization of the pump laser pulses and are induced by non-equilibrium spin-orbit 
coupled photo-holes. 

In current induced STT, spin-polarized carriers are electrically injected into a magnetic 
object, such as thin ferromagnetic layer or domain wall, from another part of a non-uniform 
magnetic structure. The physical origin of STT is in the angular momentum transfer from 
the injected carrier spins to the magnetic moments. The current induced SOT, on the other 
hand, is observed in uniform magnets with no external source of spin polarized carriers. The 
non-equilibrium spin polarization of carriers producing SOT results from current induced 
redistribution of carrier states in the band structure of the magnet. The physical origin of 
SOT is the spin-orbit coupling in the carrier bands. While the seminal works on current 
induced STT are more than 15 years ol d 12 ' 13 1 and the effect already plays a key role in 
commercially developed spintronic technologies, the research of the relativistic SOT is still at 
its infancy. Yet, the remarkable property of this inverse magneto-transport effect, allowing a 
single piece of magnet to excite itself under applied electric field, has already found practical 
utility. For example, when combined with the self-detection of the magnetization variations 
by anisotropic magnetoresistance, which is a direct magneto-transport effect based also on 
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spin-orbit coupling, an all-electric ferromagnetic resonance measurement of micromagnetic 
parameters can be performed on a single ferromagnetic nanostructure.^ 

The aim of our works reported in Ref. [3] and in this paper is to find the light induced 
torque counterparts of STT and SOT. Ferromagnetic semiconductor (Ga,Mn)As utilized in 
our experiments is a favorable candidate material for observing the phenomena. The direct- 
gap GaAs host allows the generation of high density non-equilibrium photo-carriers and the 
carrier spins interact with ferromagnetic moments on Mn via strong exchange coupling.^ 
When the ferromagnetic Mn moments are excited, this can be sensitively detected by probe 
laser pulses due to large magneto-optical (MO) signals in (Ga,Mn)As. Several groups have 
reported MO studies of fast laser induced magnetization dynamics in (Ga,Mn)As (see Sup- 
plementary information) G2H2I jj owe ver, the direct search and observation of the light induced 
torque counterparts of STT and SOT have not been the subject of these studies until our 
work in Ref. |3]and the work presented in this paper. 

In the optical spin transfer torque (OSTT), observed in our previous experiments in 
Ref. El the external source for injecting spin polarized carriers is provided by circularly po- 
larized light at normal incidence which yields high degree of out-of-plane spin-polarization 
of injected photo-carriers due to the optical selection rules in GaAs. For the photo-electrons 
producing the OSTT, the spin-lifetime is much larger than the precession time in the ex- 
change field of the Mn moments. This results in the steady state carrier spin density compo- 
nent oriented in the plane of the ferromagnetic film and perpendicular to the magnetization 
vector, which is an analogous mechanism to the adiabatic current-induced STTP Since large 
OSTT requires large spin lifetime of injected carriers, i.e. spin-orbit coupling is detrimental 
for OSTT, the weakly spin-orbit coupled photo-electrons play the key role in this caseP^ 

In the optical spin-orbit torque (OSOT), the absence of an external polarizer corresponds 
to photo-carrier excitations which are independent of the polarization of the pump laser 
pulses. Since the phenomenon relies on spin-orbit coupling, the non-equilibrium photo-holes 
are essential in this case. Because of the lack of external polarizer and because of strong spin- 
orbit coupling, precession of carriers in the exchange field of magnetic moments is not the 
origin of the torque in this case, in close analogy to the current induced SOT and reminiscent 
of the non- adiabatic STT.™ 3 The physical picture of the OSOT is as follows: The optically 
injected photo-holes relax towards the hole Fermi energy of the p-type (Ga,Mn) As on a short 
(~ 100 fs) timescale, creating a non-equilibrium excess hole density in the valence band. The 
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change in the occupation of the hole states, as compared to the equilibrium state in dark, 
and the spin-orbit coupling can generate a non-equilibrium hole spin polarization which is 
misaligned with the equilibrium orientation of Mn moments. This non-equilibrium photo- 
hole polarization persists over the timescale of the hole recombination (~ps) during which 
it exerts a torque on the Mn local moments via the kinetic-exchange coupling. The three 
key characteristic features, i.e., the non-equilibrium occupation of carrier states, spin-orbit 
coupling in carrier bands, and the resulting non-equilibrium carrier polarization acting on 
the magnetic moments via exchange coupling makes OSOT the direct optical counterpart 
of the current induced SOT. Since the applied electrical drift (and relaxation) yields a 
non-equilibrium carrier occupation in the form of asymmetric redistribution on the Fermi 
surface, the current induced SOT is observed in systems with broken inversion symmetry in 
the crystal. OSOT is caused by optical generation (and relaxation) of photo-carriers without 
an applied drift and, therefore, the broken inversion symmetry in the crystal is not required. 
It is replaced by the broken time reversal symmetry, i.e., by a non-zero spin polarization of 
carrier bands in equilibrium. This is another reason why holes in the strongly exchange-split 
valence band in (Ga,Mn)As govern OSOT. 

The schematic illustration and our experimental observation of OSOT are shown in Fig. 1. 
Apart from the physically more intriguing nature of relativistic OSOT, as compared to 
OSTT, its experimental identification is complicated by the presence of thermal excitation 
mechanisms of magnetization dynamics in the case of pump-polarization independent signals 
(see Refs. 25,26 and Supplementary information). The absorption of the pump laser pulse 
leads to photo-injection of electron-hole pairs. The non-radiative recombination of photo- 
electrons produces a transient increase of lattice temperature which builds up on the time 
scale of ~ 10 ps and persists over ^ 1000 ps. This results in a quasi-equilibrium easy- 
axis (EA) orientation which is tilted from equilibrium EA. Consequently, Mn moments in 
(Ga,Mn)As start to precess around quasi-equilibrium EA, as illustrated in Fig. la, with a 
typical precession time of ~ 100 ps given by the magnetic anisotropy fields in (Ga,Mn)As. 
As discussed in detail below, EA stays in-plane and the sense of rotation within the plane 
of the (Ga,Mn)As film with increasing temperature is uniquely defined. In the notation 
introduced in Fig. lc, the change of the in-plane angle Sep of magnetization during the 
thermally excited precession can be only positive. OSOT, illustrated in Fig. lb, acts during 
the laser pulse (with a duration of 200 fs) and fades away within the hole recombination time 
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(~ ps). It causes an impulse tilt of the magnetization which allows us to clearly distinguish 
OSOT from the considerably slower thermal excitation mechanism. Moreover, the initial 
OSOT induced tilt of magnetization can yield precession angles that are inaccessible in the 
thermally induced magnetization dynamics. This provides another evidence for OSOT. 

Examples of the direct observation of the thermally governed excitation of magnetization 
at a lower pump pulse intensity of 6/0 and of the OSOT governed excitation at a higher 
intensity of 12/ are shown in Fig. Id, where Iq — 7 /xJcm 2 . We point out that these dy- 
namical MO signals are completely independent of the polarization of pump pulses, i.e., 
they are the same for any orientation of the polarization plane of linearly polarized pump 
laser pulses and they also correspond to the polarization-independent part of the MO signals 
extracted by summing the signals induced by a + and a~ circularly polarized pump pulses 
(see Supplementary information). The distinct features of OSOT, described in the previous 
paragraph, are clearly visible when comparing the two measured trajectories of magnetiza- 
tion angles. This key demonstration has been enabled by the technique which we developed 
in Ref . [271 and which translates the measured dynamical MO signals in our pump-and-probe 
experiments to the time-dependent magnetization vector trajectory. This is done without 
assuming any theoretical model for the magnetization dynamics and without using any fit- 
ting parameter. Our experimental method utilizes different dependences of the polar Kerr 
effect (PKE) and magnetic linear dichroism (MLD) on the orientation of linear polarization 
of the probe laser pulses to disentangle the contributions to the MO signal from the out-of- 
plane and in-plane components of the magnetization motion, respectively. The magnitudes 
of PKE and MLD coefficients in a particular (Ga,Mn)As sample are determined from static 
MO experiments in which an external magnetic field is used to align the magnetization in a 
defined orientation (see Supplementary information). To obtain each point on the trajectory 
of the magnetization vector excited in the pump-and-probe experiment, we performed a set 
of measurements of the MO signal as a function of the orientation of the polarization plane 
of the linearly polarized probe pulse. The dynamical MO measurements were performed at 
zero magnetic field; prior to the experiment, the magnetization was aligned with EA. (For 
more details on the experimental technique see Refs. I3ll27l and Supplementary information.) 

Our observation of OSOT stems from an extensive growth, characterization, and MO 
study of a series of (Ga,Mn)As/GaAs materials. The ferromagnetic semiconductor epilayers 
have individually optimized molecular-beam-epitaxy and post-growth annealing procedures 
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FIG. 1: Schematic illustration and experimental observation of the optical spin-orbit 
torque, a, Schematic illustration of the thermally excited precession of magnetization M(t) around 
the transient quasi-equilibrium easy axis (EA). Mq is the magnetization vector aligned with in- 
plane equilibrium EA before the pump pulse, b, Schematic illustration of the OSOT induced by the 
in-plane transverse component s^ of the non-equilibrium hole spin polarization. On the time-scale 
of magnetization precession, OSOT causes an instantaneous tilt of the magnetization M(ti) which 
allows us to clearly distinguish OSOT from the considerably slower thermal excitation mechanism. 
The initial OSOT induced tilt of magnetization can yield precession angles that are inaccessible in 
the thermally induced magnetization dynamics, c, Definition of the coordinate system, d, Time 
evolution of the magnetization vector measured in a (Ga,Mn)As material with nominal Mn-doping 
x — 3%. The direction of the time increase is depicted by arrows. Magnetization tilt angles 5cp 
and 59 are measured with respect to equilibrium EA. Sample base temperature before pump pulse 
was 15 K and experiments were performed at zero magnetic field. At lower pump intensity 6 To 
(Jo = 7 /iJcm -2 ) the precession is induced thermally while at 12/q the OSOT induced initial tilt 
is observed. 

for each nominal Mn doping in order to minimize the density of compensating defects and 
other unintentional impurities and to achieve high uniformity of the epilayers. Nominal 
Mn-dopings in this set of 20 nm thick ferromagnetic (Ga,Mn)As epilayers span the whole 
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range from 1.5% up to ~ 14% (i.e. up to ~ 8% of uncompensated Mnc a impurities) with 
ferromagnetic transition temperatures reaching 188 K. All samples within the series have re- 
producible characteristics with the overall trend of increasing Curie temperature, increasing 
hole concentration, and increasing magnetic moment density with increasing nominal Mn 
dopingP3 The samples are in-plane magnets in which the biaxial anisotropy, reflecting the 
cubic symmetry of the host crystal, competes with an additional uniaxial anisotropy. The 
biaxial anisotropy dominates at very low dopings and the easy axis aligns with the main 
crystal axis [100] or [010]. At intermediate dopings, the uniaxial anisotropy is still weaker 
but comparable in magnitude to the biaxial anisotropy. In these samples the two equilibrium 
easy- axes are tilted towards the [110] direction and are sensitive to small changes in doping 
or temperature. At very high dopings, the uniaxial anisotropy dominates and the system 
has one strong easy- axis along the [1 10] in-plane diagonal. 

In Fig. 2 we show examples of measured dynamical MO signals in several samples with 
different magnetic anisotropics. In agreement with the doping trends in static magnetic 
anisotropics described above, we do not observe laser induced precession of magnetization 
in the very low and very high doped samples in which the magnetic easy axis aligns with one 
of the high symmetry crystal direction and its direction is insensitive to small perturbations. 
On the other hand, we observe precessions in samples with intermediate doping in which the 
EA direction is sensitive to the sample temperature and the hole concentration. Note that 
the observed precession frequencies in the studied set of samples correspond to magneto- 
crystalline anisotropy fields which are fully consistent with the respective anisotropy fields 
obtained from magnetization measurements by the superconducting quantum interference 
device (SQUID). 

The MO signals which depend on the helicity of circularly polarized pump laser pulses, 
labeled in Fig. 2 as a + — a - , originate from the excitation of ferromagnetic Mn moments 
by OSTT. The key signatures of the OSTT induced MO signals are extensively discussed 
in Ref. [3l The aim of the present work is the observation of OSOT and we therefore focus 
on magnetization precession signals which do not depend on the polarization of pump laser 
pulses. The corresponding data are labeled as a + + a~ in Fig. 2. Since the misalignment 
of the non-equilibrium hole polarization with the equilibrium orientation of Mn moments 
has the same physical origin as the dependence of the static magnetic easy-axis orientation 
in (Ga,Mn)As materials on hole density, the observation of OSOT requires a (Ga,Mn)As 
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FIG. 2: Laser pulse-induced precession of magnetization in a series of optimized 
(Ga,Mn)As epilayers. a-f, Dynamical MO signals in samples with nominal Mn-doping x — 
1.5 — 14% induced by circularly (a + and a~) polarized pump pulses were measured from which 
the helicity-sensitive [(cr + — a~)/2] and polarization-insensitive + cr~)/2] parts of the signals 
were extracted. Pump intensity in these measurements was 4/q, sample base temperature before 
pump pulse was 15 K, and external magnetic field of 20 mT was applied along the [010] crystal 
axis. Insets illustrate variations of easy axis positions across the series of measured (Ga,Mn)As 
samples at 15 K. 

film which shows sensitivity of the easy-axis orientation to small changes in doping.^ For 
measurements shown in Fig. Id and below in Figs. 3,4 we have singled out a material from 
the lower doping end in our series. This 3% Mn doped epilayer (with a Curie temperature of 
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77 K) is still a relatively low hole- density material but with already competing biaxial and 
uniaxial anisotropies for which we can expect sizable OSOT at accessible laser intensities. 

Detailed measurements of the pump-polarization independent magnetization trajectories 
for several intensities of the pump pulses are shown in Fig. 3. The key characteristics 
of the magnetization dynamics at low intensities Iq and 6/0 reflect the thermal excitation 
mechanism described in Fig. la. In equilibrium, EA in the 3% Mn doped sample is tilted by 
approximately 10° from the [010] (cp = 90°) crystal axis towards the [110] (cp = 135°) in-plane 
diagonal direction. With increasing temperature, the easy-axis rotates further towards the 
[1 10] direction. This is because the uniaxial anisotropy component scales with magnetization 
as ~ M 2 while the biaxial component scales as ~ M 4 and, therefore, the uniaxial anisotropy 
gets enhanced relative to the biaxial anisotropy with increasing temperatures. This expected 
EA rotation is confirmed by independent SQUID measurements and microscopic calculations 
based on the k • p kinetic-exchange Hamiltonianj 14 1 29 1 shown in Fig. 4a, b. The amplitude of 
the precession angles for the intensity 6/0 (Fig. 3b) is larger than for the intensity Iq (Fig. 3a). 
This is consistent with a larger increase of the transient temperature (and corresponding 
larger tilt of quasi-equilibrium EA) for the larger pump pulse intensity. We have deduced the 
temperature increase due to pump pulses from the measured precession frequencies which 
reflect the temperature dependent magnetocrystalline anisotropy energies. In Fig. 4c we 
plot the dependence of the precession frequency on the base sample temperature at low 
excitation intensity Iq and on the laser intensity at low base temperature of 15 K. From the 
comparison of these two measurements we infer the magnitude of the transient temperature 
change as a function of the laser intensity. (Note that consistent temperature vs. intensity 
calibration is obtained from the comparison of the intensity dependence of the pump-induced 
demagnetization and the temperature dependence of the remanent magnetization measured 
by SQUID.) Fig. 4c confirms a sizable difference in transient temperatures for intensities Iq 
and 6/0 • 

Remarkably, the heating of the sample by the laser pulses saturates at approximately 
10/o, as seen from Fig. 4c. The measured trajectories of the dynamical magnetization vector, 
however, show dramatic differences below and above WIq. The impulse tilt and precession 
angles inaccessible by the thermal excitations, seen in Figs. 3c, d for intensities 12/ and 30/o, 
were already pointed out in Fig. 1 as key signatures of OSOT. The complete saturation of the 
transient temperature increase at 10/q provides another confirmation that the magnetization 
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FIG. 3: Direct experimental reconstruction of the magnetization trajectory from 
MO signals, a-d, Time evolution of the magnetization vector measured in the 3% Mn-doped 
(Ga,Mn)As at pump intensities Iq — 30Iq. The plotted data correspond to time delays between 
pump and probe pulses in steps of 2 ps. Sample base temperature before pump pulse was 15 K 
and experiments were performed at zero magnetic field. 

dynamics at high pump pulse intensities are governed by a distinct non-thermal mechanism. 
The connection between OSOT and photo-carrier generation is evidenced in Fig. 4d. Here 
we show the observed change in the measured reflectivity of the (Ga,Mn)As film which 
correlates with the number of generated photo-carriers (see Supplementary information). 
The pump-induced change of the index of refraction is linear in pump-pulse intensity up to 
~ 25/o after which it starts to saturate. It means that, unlike the transient temperature, the 
number of generated photo-carriers keeps increasing with increasing pump pulse intensity 
above WIq. The concentration of photo-injected carriers can be estimated from the laser 
spot size, photon energy, and absorption and reflection coefficients of the sample. For 
the higher intensity measurements, the obtained photo-carrier density is of the order of 
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FIG. 4: Characterization of the 3% Mn-doped (Ga,Mn)As. a, Temperature dependence 
of the equilibrium EA orientation determined from SQUID magnetization measurements. Consis- 
tent EA orientations are inferred from the MO experiments, b, Microscopic calculations of the 
temperature and hole density dependent EA orientation, c, Frequency of precessing Mn moments 
measured at zero magnetic field as a function of the base temperature at low excitation intensity 
To, and as a function of the pump intensity at low base temperature of 15 K. d, Dynamics of 
pump-induced reflectivity change at intensity QIq. Inset shows the intensity dependence of the 
initial reflectivity change at base temperature of 15 K (line depicts a linear dependence). 

In Fig. 4b we show that equilibrium EA is sensitive to the hole density variations and that 
the sense of the tilt of EA with increasing hole density can be opposite than in the case of 
the temperature increase. Although not explaining the experiments on a quantitative level, 
Fig. 4b provides the clue why the instantaneous out-of-plane component of the tilt of the 
magnetization due to OSOT at higher pump intensities is opposite than the initial out-of- 
plane component of the precessing magnetization around thermally excited quasi-equilibrium 
EA at lower pump intensities. A full quantitative theory of OSOT is a challenging problem 
compared to the theory of OSTTP^ j n the latter case, the non-equilibrium spin-density of 
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weakly spin-orbit coupled photo-electrons, producing the OSTT, is directly determined by 
the external polarizer, i.e., by the intensity, propagation axis, and helicity of the circularly 
polarized pump laser beam. The relation between the non-equilibrium density of photo- 
holes and the transverse component of their spin polarization, producing OSOT in the case 
of light-polarization independent excitations, results from a more complex interplay of spin- 
orbit coupling and photo-generation and relaxation processes in the spin-split valence band 
of the ferromagnetic semiconductor. The effective field Js generating OSOT, 

(1M J~ 

— = -Mxs, 1 
dt h w 

where s is the non-equilibrium hole spin polarization density and J ~ 50 meVnm 3 is the 
hole - Mn moment exchange coupling constant, can be related in a simplified picture to the 
hole-density dependent magnetocrystalline anisotropy field (see Supplementary information 
for more theory details). Using again the k • p kinetic-exchange Hamiltonian H = Hh ost + 
JN M nS~M • cr, 1 14 , 29 1 where Uhost is the host semiconductor Hamiltonian, N Mn is the Mn local 
moment density, S = 5/2 is the local moment spin, M is the local moment unit vector, and 
a is the hole spin operator, the anisotropy field can be written as, 8 

1 d 



H, 



an 



N Mn SdM 



V f dke aM f a , k = - 1 V ( dk(a,k|-^|a,k)/ a)l 
J ^ Mn b J oJVL 



J dk(a,k\Ja\a,k)f aM = -Js. (2) 



Here e a ^ and f a ^ are the eigenenergy of H and Fermi distribution function, respectively, 
labeled by the band and wavevector index. Consistent with the in-plane orientation of 
EA, we obtain that the out-of-plane transverse component of the anisotropy field, H an ,o = 
— Jsq = for any in-plane M and any considered hole density. The in-plane transverse 
component, H arw = — Js^, is zero when M is aligned with EA at a given hole density and 
non-zero for other orientations of M at the same hole density. Since the EA orientation is 
sensitive to the hole density, as shown in Fig. 4b, H a7ljV? for a given orientation of M can 
change when the hole density is increased by, e.g., the photo-excitation (see Supplementary 
information). The sign of the calculated H an ^ is consistent with the sense of the initial 
magnetization tilt observed in experiments governed by the OSOT. The amplitudes of Js^, 
obtained from the calculated hole density dependent anisotropy fields, are ~ /xeV which is 
about 10 times smaller than the experimental strength of OSOT fields inferred form the 
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measured out-of-plane tilts of the magnetization. The k • p kinetic-exchange model is known 
to underestimate by an order of magnitude anisotropy fields in the lower doped ferromagnetic 
(Ga,Mn)As samples. 29 Considering this general limitation of the model Hamiltonian theory 
we can conclude that the calculations confirm our experimental observation of OSOT. 

To conclude, following the experimental discovery of OSTT induced by circularly polar- 
ized laser pulses, we have reported in the present paper the experimental observation of 
OSOT. The presence of this intriguing phenomenon, in which non-equlibrium photo-carrier 
polarization and the corresponding torque on magnetization are generated by spin-orbit 
coupling in the absence of an external polarizer, is evidenced by our direct experimental 
measurement of magnetization vector trajectories. OSOT which induces magnetization ex- 
citations independent of the polarization of pump laser pulses is observed in our experiments 
on picosecond timescales which is orders of magnitude shorter than in the current induced 
spin-orbit torque experiments. 
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EXPERIMENTS PRECEEDING OUR OBSERVATIONS OF OSTT AND OSOT 

(Ga,Mn)As and related ferromagnetic semiconductors are potentially ideal test-bed 
materials for exploring laser induced excitations of magnetization. Their direct-gap band 
structure allows for strong optical excitations of the electronic system, the photo -carriers can 
directly interact with magnetic moments via strong exchange coupling, the carrier mediated 
ferromagnetism produces large and tunable magnetic and magneto-optical effects, and the 
relatively simple band structure is favorable for identifying microscopic physical origins of 
the phenomena. Prior to our observations of optical spin transfer torque (OSTT) and optical 
spin-orbit torque (OSOT) presented in Ref. rllb and in this work, femtosecond laser pulse 
induced precession of magnetization in ferromagnetic (Ga,Mn)As has been reported by 
several groups [rl-r8]. In Ref. r3, magnetization precession was induced by laser pulses of a 
relatively weak intensity « 1 - 10 juJcm and ascribed to heating effects. On the other hand, 
in Refs. rl and r7 the observed precession triggered by comparably weak laser intensities was 
attributed to an effect of photo-injected holes on the magnetic anisotropy. Apart from these 
competing views, the photo-carrier based interpretation cannot be reconciled with the 
theoretical and experimental understanding of static magneto-crystalline anisotropics in 
(Ga,Mn)As [r7]. The interpretation would imply sizable changes of magnetic anisotropy in 
materials with equilibrium hole densities « 10 20 -10 21 cm" 3 induced by photo-injected holes 
of density as low as « 10 16 cm 3 . To compare with, e.g., electrical gating experiments [r9, 
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rlO], we recall that changes in the magnetic anisotropy are detected for field-induced hole 
accumulation or depletion of at least « 10 18 - 10 19 cm" 3 . We also point out that the original 
interpretation in Ref. rl in terms of photo-carrier induced changes of in-plane anisotropy was 
subsequently revised by assuming [r7,r8] out-of-plane tilts of the easy-axis. Since out-of- 
plane easy-axis rotation is not consistent with measured trends in static magnetic anisotropics 
of the considered (Ga,Mn)As/GaAs materials this again illustrates that the non-thermal laser 
induced magnetization precession due to optically generated carriers in ferromagnetic 
semiconductors has remained in Refs. rl-r8 an attractive yet unproven concept. In Ref. lib 
and in this work we provide both the direct experimental evidence for the non-thermal effects 
and the microscopic interpretation of the physical origin of the observed optical torques. 



SAMPLES 

The time-resolved magneto-optical experiments described in the main text were 
performed in a large set of optimized (Ga,Mn)As epilayers. Ferromagnetic (III,Mn)V 
semiconductors, with (Ga,Mn)As as the most thoroughly investigated example, are materials 
grown under highly non-equilibrium conditions due to the extreme levels of Mn-doping. 
Therefore, a special care has to be taken when preparing the samples and when generalizing 
the experimental results obtained in a given sample to the universal behavior of this material 
system. Very recently we reported a systematic study of optical and magneto-optical 
properties of an optimized set of (Ga,Mn)As epilayers spanning the wide range of accessible 
substitutional Mno a dopings [rl 1]. The optimization of the materials in the series, which is 
performed individually for each nominal Mn doping, minimizes the uncertainties in the 
experimental sample parameters and produces high quality epilayers which are as close as 
possible to uniform uncompensated (Ga,Mn)As mixed crystals. For each nominal Mn doping 
x, the growth and post-growth annealing conditions were separately optimized in order to 
achieve the highest Curie temperature T c attainable at the particular x. The highest T c criterion 
was found to lead simultaneously to layers with maximized uniformity and minimized 
compensation by unintentional impurities and defects [rl 1]. All samples are in-plane magnets 
in which the cubic anisotropy competes with an additional uniaxial anisotropy. At very low 
dopings, the cubic anisotropy dominates and the easy axis aligns with the main crystal axis 
[100] or [010]. At intermediate dopings, the uniaxial anisotropy is still weaker but comparable 
in magnitude to the cubic anisotropy. At very high dopings, the uniaxial anisotropy dominates 
and the system has one strong easy-axis along the [1-10] in-plane diagonal. We confirmed 
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that the laser pulse-induced precession of magnetization can be observed in all (Ga,Mn)As 
epilayers except in those with very low and very high doping levels where one of the 
anisotropics strongly dominates - see Table I. 



sample 
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11.5 


Yes 


F002 


3 


77 


16.2 


Yes 


F016 


3.8 


96 


24.7 


Yes 


E101 


4.5 


111 


27.8 


Yes 


F020 


5.2 


132 


33.3 


Yes 


E115 


7 
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179 


63.7 
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F056 


14 


182 


78.1 


No 



TABLE I: Table summarizing the basic characteristics of selected samples from the series of optimized 
materials: x is nominal doping, T c is Curie temperature, M s is saturated magnetic moment. Last column shows if 
the impact of ultrafast laser pulse induces a precession of magnetization for the sample temperature of 15 K and 
the pump pulse fluence of about 30 jiJ.cm" 2 . 

For simplicity, we show in this Supplementary information only the helicity- 
independent signals measured in (Ga,Mn)As epilayer with a nominal doping x = 3%, which is 
the sample that was used for a detailed analysis presented in the main text. The extensive 
information about the helicity-dependent signals can be found in our previous paper [rl lb]. 

EXPERIMENTAL DETAILS 

We investigated the laser-pulse induced dynamics of magnetization by a pump-and- 
probe magneto-optical technique. A schematic diagram of the experimental set-up is shown 
in Fig. SI. The output of a femtosecond laser is divided into a strong pump pulse and a weak 
probe pulse that are focused to the same spot on the measured sample. Laser pulses, with the 
time width of 200 fs and the repetition rate of 82 MHz, were tuned to 1.64 eV, i.e. above the 
semiconductor band gap, in order to excite magnetization dynamics by photon absorption. 
The pump pulses were usually circularly polarized (with a helicity controlled by a wave plate) 
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and the probe pulses were linearly polarized. The measured magneto-optical signals 
correspond to the probe polarization rotation induced by the pump pulses (see Fig. SI). The 
experiment was performed close to the normal incidence geometry = 2° and 8° for pump 
and probe pulses, respectively) with a sample placed in a cryostat, which was placed between 
the poles of an electromagnet. The external magnetic field H ext was applied in the sample 
plane at an angle <^#with respect to the [100] crystallographic direction in the sample plane 
(see Fig. SI). 




Fig. S 1 . Schematic diagram of the experimental set-up for a detection of the magnetization precession induced 
in (Ga,Mn)As by an impact of the circularly polarized femtosecond laser pump pulse. Rotation of the 
polarization plane of reflected linearly polarized probe pulses is measured as a function of the time delay A* 
between pump and probe pulses. The orientation of magnetization in the sample is described by the polar angle cp 
and azimuthal angle 6. The external magnetic field H ext is applied in the sample plane at an angle <p H . 

DETERMINATION OF THREE-DIMENSIONAL MAGNETIZATION VECTOR 
TRAJECTORY FROM DYNAMICAL MAGNETO-OPTICAL SIGNALS 

In (Ga,Mn)As there are two magneto-optical (MO) effects that are responsible for the 
measured signal [rl2]. In the following we will concentrate on the rotation of the polarization 
plane A/3 of the reflected linearly polarized light (see Fig. S3(b) for a definition of the probe 
pulses polarization plane orientation /J) but the same apply also for the change of the light 
ellipticity (see below). We will also limit the discussion to the case when the incident light is 
close to the normal incidence. The first of the MO effects is the well-known polar Kerr effect 
(PKE), which is sometimes called magneto-optical Kerr effect (MOKE), where the rotation of 
polarization occurs due to the different index of refraction for a + and a ' circularly polarized 
light propagating parallel to the direction of magnetization - see Fig. S2(a). Consequently, Af3 
due to PKE is proportional to the projection of magnetization to the direction of light 
propagation. We also note that this MO effect is linear in magnetization (i.e., the sign of the 
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polarization rotation is changed when the direction of magnetization is reversed) and that A/5 
is independent on (5 [r!2]. 




Fig. S2. Schematic illustration of two magneto -optical effects that are responsible for a rotation of the 
polarization plane Aft of reflected light at normal incidence in (Ga,Mn)As. (a) Polar Kerr effect (PKE) that is due 
to the different index of refraction for cr + and <j~ circularly polarized light propagating parallel to the direction of 
magnetization M. (b) Magnetic linear dichroism (MLD) that is due to the different absorption (reflection) 
coefficient for light linearly polarized parallel and perpendicular to M if the light propagates perpendicular to the 
direction of M [rl2]. 

The second MO effect is the magnetic linear dichroism (MLD), which originates from 
different absorption (reflection) coefficient for light linearly polarized parallel and 
perpendicular to M, that occurs if the light propagates perpendicular to the direction of 
magnetization M - see Fig. S2(b) [rl2]. Consequently, A/5 due to MLD is proportional to the 
projection of magnetization to the direction perpendicular to the direction of light 
propagation, A/3 is quadratic in magnetization (i.e., the sign of A/3 is not changed when the 
direction of magnetization is reversed) and it varies as sin(2/?) [rl2]. 

For normal incidence of light, the simultaneous presence of two MO effects that are 
sensitive to the out-of-plane orientation of magnetization (PKE) and the in-plane orientation 
of magnetization (MLD) enables to perform a full three-dimensional stroboscopic 
reconstruction of magnetization vector movement from the measured dynamical MO signals 
[rl2]. In particular, different polarization dependences of PKE and MLD enable us to 
disentangle the pump-induced out-of-plane and in-plane motions of magnetization, 
respectively [rl2]. In equilibrium, the magnetization points to the easy axis direction, which 
can be characterized by the polar angle <po and azimuthal angle do (see Fig. S3(b) for a 
definition of the coordinate system). The impact of the pump laser pulse induces a 
precessional motion of magnetization. The magnetization orientation at any time delay At is 
given by (p(Ai) = cpo + S(p(At) and 0(At) = 9 + S9(At). The time-dependent position of 
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magnetization leads to a time-dependent MO response of the sample that, moreover, depends 
on the orientation of polarization plane of probe pulses /? - i.e., the measured dynamical MO 
signals depend both on At and /?: 5MO{At, p). Consequently, it is possible to determine a time 
evolution of the magnetization vector orientation - without any numerical modeling - directly 
from the measured MO data [rl2]. In particular, 

Sp(At) = [5MO(At,<p )-SMO(At,<p -90 o )]l{4P MLD ), (1) 

S0(At) = -[5MO(At,(p -45°)+mo(At,<p -135°)]/(2P P ^). (2) 

where 5MO{At, (po) is the pump-induced dynamical MO signal measured by probe pulses with 
P = (po\ 5MO{At, cpo - 90°) is the pump-induced dynamical MO signal measured by probe 
pulses with (3 = <p - 90°; etc. P PKE and P MLD are MO coefficients that describe the magnitude 
ofPKEand MLD. 

We note that for a quantitative evaluation of the magnetization motion from the 
measured MO signals it is essential to know the magnitude and the sign of the corresponding 
MO coefficients P PKE and P MLD in the investigated sample [rl2]. These coefficients can be 
directly measured if the magnetization is oriented by a strong external magnetic field to the 
requested position. For a measurement of P PKE , the magnetization has to be oriented 
out-of-plane (i.e., 6- 0°) and any orientation of light polarization plane /? can be used. For a 
measurement of p MLD ? the magnetization has to be oriented in-plane (i.e., = 90°) and P = 

PKE rdMLD 

cpo - 45° has to be used [rl2]. In Fig. S3(a) we show the spectral dependence of P and P 
measured in a (Ga,Mn)As epilayer with nominal doping x = 3% that is the sample used for a 
detailed study of the magnetization dynamics at different excitation intensities, as described in 
the main paper. 
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Fig. S3, (a) Spectral dependence of static PKE and MLD in (Ga,Mn)As epilayer with nominal doping x = 3%, 
the arrow indicate the spectral position of laser pulses used in the time -resolved experiment, (b) The orientation 
of magnetization M in the sample is described by the polar angle q> and azimuthal angle the orientation of the 
polarization plane E of probe pulses is described by the angle /?. 

In Fig. S4 we show MO hysteresis loops that were measured for the in-plane 
orientation of the external magnetic field (see the inset in Fig. S4). The observed M-shaped 
hysteresis loop is a signature of the existence of four energetically equivalent easy axis 
positions of magnetization in a epilayer. Prior to the time-resolved experiment, we prepared 
the magnetization in a state close to the [010] crystallographic direction (i.e., in an easy axis 
position labeled "1" in Fig. S4). 
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Fig. S4. Magneto-optical hysteresis loops measured in a (Ga,Mn)As epilayer with nominal doping x = 3 % at 
15 K; the four energetically equivalent easy axis positions of magnetization are schematically labeled "1" to "4" 
in the inset. 

The impact of a laser pulse on (Ga,Mn)As induces a precession of magnetization by 
two distinct mechanisms - the helicity-dependent and the helicity-independent [rllb]. To 
separate these effects, we computed from the dynamical MO signals measured after excitation 
by cr + and cr~ circularly polarized pump pulses the helicity-independent [(cr + + cr~ )/2] and 



helicity-dependent [(cr + - a~ )/2] signals. In Fig. S5(a) we show the as-measured data for cr + 
and a~ excitation. The helicity-independent part of the signal, which is the main focus of this 
paper, is shown in Fig. S5(b). We note that the helicity-independent signal is identical to that 
measured with linearly polarized pump pulses [with any orientation of the polarization plane - 
see Fig. S5(b)]. 



30 



^ 20 
1 10 



i 




(a) 



3 



500 1000 1500 2000 2500 
Time delay (ps) 



a + a 

s linear polarization 
p linear polarization 




500 1000 1500 2000 2500 
Time delay (ps) 



Fig. S5. (a) Dynamical magneto-optical (MO) data measured for circularly (cr + , a ) polarized pump pulses, (b) 
Comparison of the helicity-insensitive [(a + + a ~ )/2] part of the signals shown in (a) with the MO data 
measured for linearly polarized (s, p) pump pulses. Excitation intensity I = 7 jaJ.cm 2 , sample temperature T = 
15 K, probe polarization orientation J3 = 0°, no external magnetic field applied. 

The dependence of helicity-independent dynamical MO signal on the probe 
polarization orientation /? is shown in Fig. S6(a). The measured data can be fitted well by the 
phenomenological equation 



SMO(t) = Acos(cd m J + A)e / T(J + Ce ' Tp , 



(3) 



where A and C are the amplitudes of the oscillatory and pulse function, respectively, co Mn is 
the ferromagnetic moment precession frequency, A is the phase factor, x G is the Gilbert 
damping time, and z p is the pulse function decay time. All the measured data in Fig. S6(a) can 
be fitted well by Eq. (3) with a one set of parameters cdmu, ?g and r p . The dependence A{(3) 
obtained by this fitting procedure is displayed in Fig. S6(b). The position of the maximum in 
the dependence A(/S) at (5 « 100° corresponds to the equilibrium position of magnetization in 
the sample cpo [rl2]. We recall that prior to this measurement we prepared the magnetization 
in the easy axis labeled "1" in Fig. S4, which close to [010] crystallographic direction. 
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Time delay (ps) p (deg) 

Fig. S6 (a) Dynamics of the helicity-independent MO signal induced by an impact of pump pulse on the sample 
that was measured by probe pulses with different probe polarization orientations p (points), (b) Probe- 
polarization dependence of the amplitude of the oscillatory part A that was obtained by fitting the measured MO 
dynamics by Eq. (3) [lines in part (a)] with a fixed values co Mn = 33.9 GHz, t g = 165 ps, and x p = 880 ps (points). 
The line is a fit (see Ref. r 12 for details). The vertical arrow depicts the derived easy axis position in the sample. 
Excitation intensity I = 7 jaJ.cm" 2 , sample temperature T= 15 K, no external magnetic field applied. 

In Fig. S7 we show the dynamics of magnetization that was deduced from the data 
depicted in Fig. S6(a) using Egs. (1) and (2) - in part (a) and (b) the time evolutions of Sep and 
SO, and the corresponding polar plot are depicted, respectively. 




Time delay (ps) 59 (deg) 



Fig. S7 Magnetization real-space trajectory deduced from the MO data shown in Fig. S6. (a) Time evolution of 
Sep and 86. (b) Orientation of magnetization at different times after the impact of the pump pulse; the equilibrium 
position of the easy axis is depicted by the grey spot. 

In principle, the measured dynamical MO signals can contain information not only 
about the pump-induced magnetization change but also about the pump-induced change of the 
complex index of refraction (the so-called "optical part" of the MO signal) [rl3]. 
Consequently, the dynamics of both the rotation and ellipticity has to be measured and 
compared before the obtained magneto-optical signal is attributed to the magnetization 
dynamics [r!3]. In Fig. S8 we show the corresponding data measured at pump intensity 18/# = 
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126 |uJ.cm" 2 , which is in the middle of the investigated intensity range (see the main paper). 
The measured dynamics of the rotation and ellipticity are identical that confirms that the 
procedure described above can be meaningfully used for a reconstruction of the magnetization 
movement. In particular, we would like to note that the displayed data which were measured 
for f5 - 105° « (po provide direct information about the in-plane, S<p, movement of the 
magnetization [cf. Eq. (1)]. At low excitation intensities, the MO data measured for this value 
of P show a monotonous rise of the (negative) dynamic MO signal for time delays < At < 50 
ps - see Fig. S6(a) for intensity Io = 7 |uJ.cm" 2 . Taking into account that F^ LD < in the used 
spectral range [see Fig. S3(a)], this leads to S(p>0 - i. e., the magnetization starts its 
movement towards the [-110] crystallographic direction in the sample (see Fig. S7). On the 
contrary, at high excitation intensities the measured MO signal has a opposite sign for a time 
delay < At < 50 ps - see Fig. S8 for intensity 18Io = 126 |uJ.cm" 2 - and, therefore, Sep < in 
the initial stages of the magnetization movements - i. e., the magnetization starts its 
movement towards the [010] crystallographic direction in the sample (see also Fig. 3 in the 
main paper). 




Time delay (ps) 



Fig. S8. Comparison of the laser pulse -induced change of the rotation and ellipticity; the data are normalized. 
Excitation intensity / = 18/ = 126 jiJ.cm" 2 , sample temperature T = 15 K, probe polarization orientation /? = 
105°, no external magnetic field applied. 
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CALIBRATION PROCEDURE USED FOR EVALUATION OF THE PUMP- 
INDUCED TEMPERATURE AND HOLE CONCENTRATION INCREASE 

Absorption of pump laser photons leads to a photo-injection of electrons and holes 
with identical concentrations. However, because the density of states is considerably larger 
for holes (due to their larger mass), the pump-induced change of the distribution function is 
much larger for electrons and, consequently, the measured change of complex index of 
refraction is dominated by electrons [rl4]. Therefore, the dynamics of the reflectivity change 
AR/R - that is shown in Fig. S9(a) - provide information about the lifetime of photo -injected 
electrons in the investigated sample. From the measured data we can conclude that the 
population of photo-injected free electrons decays with a time constant (defined as a decrease 
to 1/e) of « 15 ps. This rather short lifetime of free electrons is similar to that reported for the 
low temperature grown GaAs (LT-GaAs), which is generally interpreted as a result of a high 
concentration of nonradiative recombination centers induced by the low temperature growth 
mode of the MBE [rl5]. The nonradiative recombination of photo-injected electrons is 
accompanied by an emission of phonons (i.e., by a heating of the sample). Consequently, the 
measured decay time of AR/R corresponds to a rise time of the laser-induced transient change 
of the sample temperature ST. We note that the lifetime of the photo-injected holes can be 
expected to be longer than that of electrons because (Ga,Mn)As is a heavily p-doped material 
and, therefore, the traps for holes are at least partially saturated by the background holes. 

We applied the following analysis for an evaluation of the laser-induced change of the 
hole concentration dp. The initial value of AR/R can be regarded as a measure of the relative 
number of the photo-generated electrons (and thus also of the holes). It is apparent from the 
inset in Fig. S9(a) that the pump-induced change of the complex index of refraction is linear 
in a pump-pulse intensity up to / « 25Io =180 |uJ.cm" 2 after which the initial values of AR/R 
start to saturate. The thickness of the investigated (Ga,Mn)As film (20 nm) is considerably 
smaller than the penetration depth of the laser at the corresponding photon energy (« 600 nm). 
Consequently, the absorption of photons leads to a rather constant concentration of photo- 
injected carriers in (Ga,Mn)As along the growth direction, which can be estimated from the 
laser spot size on the sample, photon energy, and absorption and reflection coefficients of 
GaAs. On the other hand, also the carriers photo -injected in the substrate could, in principle, 
contribute to the measured signals but we have verified that the measured signals in 
(Ga,Mn)As samples on GaAs substrate are markedly different from that in a bare GaAs 
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substrate. In Fig. S10 we show the estimated intensity dependence of the initial concentration 
of the excess holes dp. 
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Fig. S9 (a) Dynamics of pump pulse-induced reflectivity change AR/R at intensity 6I = 42 jaJ.cm" ; Inset: 
Intensity dependence of the initial value of the reflectivity change at 15 K (points), line depicts the linear 
dependence, (b) Frequency of precessing Mn moments measured (without external magnetic field) as a function 
of base temperature at low excitation intensity I = 7 jaJ.cm" 2 and as a function of the laser intensity at low base 
temperature of 15 K. This figure is re -plotted from the main paper for convenience. 



As already described above, the fast nonradiative recombination of photo -injected 
electrons induces a transient pump pulse-induced increase of the lattice temperature ST within 
tens of picoseconds after the impact of the pump pulse and this temperature increase persists 
over hundreds of picoseconds. An estimate of ST can be obtained from the measured 
precession frequencies because they correspond to magneto -crystalline anisotropy fields in 
the sample [rllb], which are rather strongly temperature dependent. In Fig. S9(b) we plot the 
dependence of the precession frequency on the base sample temperature at low excitation 
intensity Iq and on the laser intensity at low base temperature of 15 K. From the comparison 
of these two measurements we infer the magnitude of the transient temperature change ST as a 
function of the laser intensity, which is shown in Fig. S10. (We note that consistent 
temperature vs. intensity calibration is obtained from the comparison of the intensity 
dependence of the pump-induced demagnetization and the temperature dependence of the 
remanent magnetization measured by SQUID.) Remarkably, the temperature increase 
saturates at a considerably lower pump intensity than the hole concentration increase does that 
enables us to identify clearly the influence of the latter in the magnetization dynamics (see the 
main paper). 
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Fig. S10 The pump-pulse induced concentration of the excess holes dp and the transient temperature increase ST 
as a function of the laser intensity. 



QUALITATIVE THEORY OF OSOT 

Absorption of pump laser pulse with a photon energy of 1.64 eV leads to photo- 
injection of charge carriers with a rather high excess energy (the low-temperature band gap 
energy of GaAs is 1.52 eV). Immediately after photoinjection, the carriers have a non-thermal 
distribution. The Fermi-Dirac distribution is established on a time scale from tens to hundreds 
of femtoseconds (mainly due to mutual collisions between the carriers) but the resulting 
characteristic carrier temperature significantly exceeds that of the lattice. The carrier 
temperature starts to cool down towards the lattice temperature (by an emission of phonons) 
on a time scale from hundreds of femtoseconds to several picoseconds and, finally, the non- 
equilibrium carriers recombine. During this process holes relax through the available states in 
the spin-split valence band of the ferromagnetic semiconductor. Consequently, a change in the 
occupation of the hole states, as compared to the equilibrium state in dark, leads - due to the 
spin-orbit coupling - to a non-equilibrium hole spin polarization which is misaligned with the 
equilibrium orientation of Mn moments. This non-equilibrium photo-hole polarization persists 
over the timescale of the hole recombination during which it exerts a torque on the Mn local 
moments via the kinetic-exchange coupling. The exact theoretical treatment of the mutual 
interplay of all the involved effects is rather complicated and, therefore, derivation of a full 
quantitative theory of OSOT is a challenging problem. (We note that in the case of OSTT the 
non-equilibrium spin-density of weakly spin-orbit coupled photo-electrons, producing the 
OSTT, is directly determined by the external polarizer - i.e., by the intensity, propagation 
axis, and helicity of the circularly polarized pump laser beam - which make the corresponding 
theory much simpler [lib]). As a first proxy to the complex theory of the non-equilibrium 
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photo-hole spin polarization, related to OSOT by Eq. 1 in the main text, we consider the hole 
density dependent magnetic anisotropy fields (see Eq. 2 in the main text). The gross 
simplification in this qualitative theory is that we consider OSOT produced by fully 
equilibrated photo-holes acting on the local moments over the recombination time. Within this 
approximation, the misalignment of the non-equilibrium hole polarization with the 
equilibrium orientation of Mn moments has the same physical origin as the dependence of the 
static magnetic easy-axis orientation on the hole density [r!6]. 
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Fig. Sll Computed in-plane angular dependence of the anisotropy energy E (a) and the anisotropy field H an>(p = 
SEIScp (b) in (Ga,Mn)As with 3% of Mn and strain - 0.01% at K. Inset: Dependence of the easy axis 
position on the concentration of holes. 



In Fig. SI 1(a) we show the in-plane angular dependence of the anisotropy energy in 
(Ga,Mn)As, which was computed using the kp kinetic-exchange Hamiltonian [rl6]. The 
position of the energy minimum (i.e., the easy axis position) depends strongly on the 
concentration of equilibrium holes [see the inset in Fig. SI 1(a)]. The femtosecond pump-pulse 
induced change of hole concentration Sp (from the sample equilibrium hole concentration at 
dark po) leads to an abrupt in-plane movement of the energy minimum and, consequently, to a 
misalignment of the energy minimum and the equilibrium orientation of magnetization M. 
The resulting anisotropy field, which is aiming at having them aligned again, has a zero out- 
of-plane component H an ,d= SEIS6 and a non-zero in-plane component H an><p = SEIScp [see Fig. 
SI 1(b)]. Consequently, the resulting quasi-instantaneous torque acting on M is oriented in the 
out-of-plane direction (see Eq. (1) and Fig. 1(b) in the main paper) which explains why we 
observed the fast out-of-plane magnetization tilt for high excitation intensities [see Fig. 3(c) 
and (d) in the main paper]. On the other hand, at low excitation intensities the pump-induced 
temperature increase is considerably slower and the resulting temperature -related movement 
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of the energy minimum (quasi-equilibrium easy axis) is in the opposite direction than in the 
case of the hole concentration-related movement [see Fig. 4(b) in the main paper]. This 
explains why we observed at low excitation intensities a much slower rise of the out-of-plane 
component of the precessing magnetization and towards the opposite direction than in the 
case of OSOT at high excitation intensity [see Fig. 3(a) and (b) in the main paper]. The 
simplified analysis therefore provides a physically transparent qualitative description of the 
experimentally observed OSOT. 
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